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We present analytical results for the strongly anisotropic random 
field Ising model, consisting of weakly interacting spin chains. We 
combine the mean-field treatment of interchain interactions with an 
analytical calculation of the average chain free energy ("chain mean- 
field" approach). The free energy is found using a mapping on a 
Brownian motion model. We calculate the order parameter and give 
expressions for the critical random magnetic field strength below which 
the ground state exhibits long range order and for the critical tem- 
perature as a function of the random magnetic field strength. In the 
limit of vanishing interchain interactions, we obtain corrections to the 
zero-temperature estimate by Imry and Ma [Phys. Rev. Lett. 35, 1399 
(1975)] of the ground state density of domain walls (spin-flips) in the 
one-dimensional random field Ising model. One of the problems to 
which our model has direct relevance is the lattice dimerization in 
disordered quasi-one-dimensional Peierls materials, such as the conju- 
gated polymer irans-polyacetylene. 
PACS numbers: 75.10.-b, 05.50. +q, 64.60. Cn 
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I. INTRODUCTION 



Owing to the wide range of physical systems that can be described by Ising-like 
models, it is believed that the random field Ising model (RFIM), i.e., the Ising model 
with a spatially random magnetic field, captures the essential physics of many disordered 
systems. Therefore, over the past two decades this model has been extensively studied, 
both theoretically and experimentally 0. The most studied realization of the RFIM is 
the dilute antiferromagnet in a uniform field. This realization was first found by Fishman 
and Aharony || in 1979 and is still under current investigation Other equally 

interesting as diverse problems to which the RFIM has been found to be applicable are, 
e.g., the Anderson-Mott transition of disordered interacting electrons and aspects of 
protein folding 0. 

Much of the theoretical interest in the RFIM has originated from the early work by 
Imry and Ma ||. In 1975, these authors noticed that in the presence of a random magnetic 
field the lowest energy state of the Ising model may have no long range order (LRO). 
They considered the energetics of creating a "wrong" domain {i.e., a domain in which all 
spins have reversed signs) of linear size R in the ordered phase. The typical domain size 
R caused by the disorder fluctuations is determined by the balance between the energy 
cost of domain walls and the possible energy gain due to the decrease of the interaction 
energy between the reversed spins and the random field. In the one-dimensional case, 
this balance reads: 



where / is the energy of one domain wall and (h 2 ) is the strength of the random magnetic 
field. This zero temperature estimate implies that in the one-dimensional random field 
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(1) 



Ising model LRO does not exist, as in an infinite system it is always possible to choose 
R large enough to make the creation of a "wrong" domain favorable. From Eq.(|l|), the 
zero temperature density of spin-flips is estimated to be 

n s = - ~ 2 
R P y 1 



The lower critical dimension above which LRO in the RFIM is possible is 2 ||10|. 

In most previous studies of the RFIM, the interactions between the Ising spins were 
assumed to be isotropic, i.e., independent of the direction. These studies focused on the 
universal critical properties of the model, such as the critical exponents and the lower 
critical dimension, which do not depend on the degree of anisotropy of the interactions. 
On the other hand, the anisotropic RFIM has several nonuniversal properties, such as 
the critical temperature, that are very interesting in relation to disordered quasi-one- 
dimensional systems. 



For example, we recently showed [IT] that disordered Peierls chains with a doubly 
degenerate ground state can under certain conditions be described by the one-dimensional 
RFIM. In this mapping, the Ising variables describe the local lattice dimerization, while 
the random field corresponds to disorder in the electron hopping amplitudes. The spin- 
flips induced by the random field then correspond to disorder-induced neutral solitons 
in Peierls chains. The density of spin-flips is related to the magnetic properties of the 
disordered Peierls systems, since the neutral solitons have spin 1/2 ||12|| . The conjugated 
polymer ircms-polyacetylene is a well-known example of a Peierls system for which this 
mapping may be used. Taking into account the weak interchain interactions in quasi- 
one-dimensional crystalline trans-polyacetylene, the thermodynamics of solitons in this 



material may be described by the anisotropic RFIM |13J. In this context, the interesting 



properties of the anisotropic RFIM are the temperature of the three-dimensional Peierls 



transition, the order parameter, and the density of spin-flips as a function of the random 
field strength and (or) temperature. These properties strongly depend on the degree of 
anisotropy of the interactions. 

In the present paper we study the <i-dimensional anisotropic RFIM (d > 2), which 
consists of regularly arranged chains of spins in a random magnetic field (Sec. [TT]). The 
interaction between spins in neighboring chains is supposed to be much smaller than 
the interaction within chains and will be treated in the mean-field approximation; the 
intrachain interactions will be treated exactly. Using the transfer matrix formalism and 
a mapping on a Brownian motion model, we find an analytical expression for the average 
free energy of the continuum version of our model (Sec. |ITT|). This key result is used in 



Sec. [IV] to obtain the order parameter, the critical temperature, and the density of spin- 
flips. In Sec. [IV A| , we present analytical results for these quantities in the limit where 
thermal creation of spin-flips dominates their creation by the random magnetic field; 
in Sec. [IV lj| we consider the opposite limit. We find the critical random field strength 
below which the ground state exhibits LRO; well below the critical strength, the density 
of spin-flips is found to be exponentially suppressed. 

In the absence of the interchain interactions, the anisotropic model reduces to the 
one-dimensional RFIM. In this limiting case, our solution is exact. It shows that the 
estimate for the density of spin-flips Eq.(^j) in a single chain is correct to lowest order in 
the random magnetic field strength and gives higher-order corrections to this expression. 
Previous attempts to obtain an exact solution for the one-dimensional case |T^] did not 
yield closed expressions for the free energy and the spin-flip density. 

We conclude this paper by giving the physical interpretation of our results in Sec. [V] 
and summarizing in Sec. [VT] . 
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II. MODEL 



We consider the anisotropic RFIM on a <i-dimensional lattice with a total number 
of N d sites. The Ising spins occupying the lattice, interact with their nearest neighbors 
and with a magnetic field that has a random value at each lattice site. The anisotropy 
results from the fact that the exchange interaction, /, along one of the lattice directions 
is assumed to be much stronger than the interaction, I±, perpendicular to this direction: 
< <C / (for definiteness, the interactions are chosen positive). The lattice is thus 
divided into N^ 1 chains (labeled a) of N spins (labeled i). The energy of our model is 
given by: 

= I L 2 L ^ 2 L^^- ( 3 ) 

(i,i'),a ~ i,(a,a') ~ i,a 

Here, a ia = ±1 is the classical spin variable describing the two possible spin directions, 
up and down, at site i of chain a. The summations over (i, i') and (a, a') are restricted 
to nearest neighbor pairs. The random magnetic field hi a is assumed to be uncorrelated 
for different lattice sites, with vanishing mean and Gaussian correlator of strength e: 

(hux) = and (h ia h jf3 ) = e5^8 a3 . (4) 

Here, (...) denotes the average over the random field realizations. Finally, we impose 
open boundary conditions on the model. This choice does not affect our results, as these 
will be derived for a large system, iV — > oo. 

For vanishing magnetic field strength e, the ground state of the spin system Eq. @ is 
ferromagnetically ordered. This is not necessarily the case when a finite random field is 
present: For a number of spins it may then be favorable to flip and align with the local 
direction of h{ a to minimize the total energy. 
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To calculate the partition function of the model we treat the interchain interactions 
in the mean-field approximation. Then the energy Eq.(^) becomes a sum of single-chain 
energies and can be written as 



£ MF [h] = N*- 1 ^ 



- (1 - (Jia i+1 ) 



B 2 I ' 
(B + hi)a % + j- + -j 



(5) 



where B is the homogeneous mean field, I± equals I± times the number of nearest- 
neighbour chains of one particular chain, and we took into account that the contributions 
of all N d ~ Y chains are on average equal. The mean field B satisfies the self-consistency 
equation: 



B = ~M<<7», 



(6) 



where ((...)) denotes both the thermal average and the average over random field re- 
alizations. The order parameter ((c)), which is the normalized difference between the 
average number of sites with spin up (iVj) and spin down (iVj, 



iVj - N 
iV T + N 



(7) 



can be found from 



]__d_ 



(F). 



Here, (F) denotes the average over the random field realization of the free energy of a 
single chain in the mean field B. For a particular random field realization: 



-0E[a,h] 



where the chain energy E [a, h] is given by the relevant part of Eq. ( 



(9) 



E[a,h] = £ 



- (1 - <Ti<r i+1 ) - (B + hi)ai 



(10) 



Using Eq.([8]), the explicit form of the self-consistency condition (^) becomes: 

B = -^< F >- (11 > 

In view of the relative weakness of the interchain interactions (I± <C I), this so-called 



chain mean-field approximation []15] is expected to be accurate. 

The disorder-averaged free energy, (F), which we calculate analytically in the next 
section, can also be used to find other relevant thermodynamic properties. In particular, 
the average density of spin-flips inside a chain as a function of temperature and random 
field strength is obtained through: 

n. = (12) 

since I is the spin-flip creation energy. In the context of the degenerate ground state 
conjugated polymers, mentioned in the Introduction, n s directly gives the density of 
solitons within chains. 



III. THE AVERAGE FREE ENERGY 

To derive the continuum version of the discrete model Eq. (|10D , we follow the standard 
path |l6||rjj by first rewriting the chain partition function in the transfer matrix formal- 
ism. The partition function is thus considered as the expectation value of an ordered 
product of operators along the chain axis: 

TV 

Z[h] = £exp{ -PE[<r, h) } = Y, i a ' I II Ti k). (13) 

{a} a,a'=±l 1=1 

Here, the final summation over a, o' accounts for all possible boundary conditions for 
the spins on the chain's ends. On the basis spanned by the two vectors | + ) and | — ) 
(corresponding to a = +1 and a — — 1, respectively), the transfer matrix % reads: 
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Ti 



,+f3(hz+B) e -/3I e +f3(hz+B) 



e -/3I e -/3(h,+B) e -/3(hi+B) 



(14) 



Ti depends on the site-index i through the random magnetic field. 

A relation between the transfer matrix Eq.(|H|) and its (euclidean) quantum Hamil- 
tonian is established by interpreting the axis of the chain as the (imaginary) time axis 
of quantum mechanics. This means that neighboring sites labeled by i are considered as 
subsequent times t in the continuum model. In this sense carries information about 
the time evolution of the system and can in fact be identified with the time evolution 
operator T of a corresponding Hamiltonian H: 



f = exp{-#} 



(15) 



In the continuum approximation the deviation of T from the identity matrix is assumed 
to be small, so that the last expression can approximately be written as 



T = 1 — if, 



(16) 



where H, obtained from Eq . (|i~4D , has the form: 



H(t) = /3(h(t) + B)& 3 + exp{-/3I}di. 



(17) 



Here, &i and are the Pauli matrices and the random magnetic field is now Gaussian 



in time: 



(h(t)) = and (h(t)h(t')) = e5(t-t'). 



18) 



In the last expression, we have chosen the lattice constant as the unit of length (time). 
The Hamiltonian Eq . (|i~7|) describes the relaxation of a spin 1/2 in a magnetic field which 
has a constant x-component and a random and time- dependent z-component. 
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It can easily be seen from Eq.(|T4]) that the validity of the expansion Eg. (|T6|) , i.e., 
the validity of the continuum description of the discrete Ising model, requires the two 
conditions 

exp{-/5/} < 1 (19) 

and 

(3(\h(t)\ + \B\) « 1 (20) 

to be fulfilled simultaneously While the first condition results in an upper limit 
for the temperature, T <C /, the latter implies a lower temperature limit: T ^> 
O^maxj-y/e, 

The continuum limit allows us to write the partition function Eq.(|I5D as a sum over 
matrix elements for the time-ordered propagation of a spin-1/2 particle in imaginary 
time: 

Z[h] = (v'\f expi f tf H(t)dt\ \a). (21) 

<T,er'=±l ^ ' 

Here, T is the time-ordering operator and the integration over the dimensionless time- 
variable t has to be performed from the initial time ti = to the final time tf = N 
corresponding to the chain length. 

It is now convenient to rotate the spinor basis over an angle tt/2 around the 02-axis, 
which yields the transformed Hamiltonian: 

H\t) = exp ji^<7 2 j H(t) expj— i^o- 2 | 

= -P{h{t) + B)a 1 + ex V {-(3I}a 3 , (22) 

and 



Z[h] = 2( + \t explj" H'(t)dt^ |+). (23) 
To make this expression more compact, we introduce the wave function of the spin, 

¥(t) = Mt)\ + ) + ^(*)|->, (24) 

which obeys the time-dependent Schrodinger-like equation, 

= (25) 

Together with the inital condition $?(t = 0) = | + ), the partition function Eq.(^) now 
simply reads: 

Z[h] = 2^(N). (26) 

Our problem of calculating the partition function has, of course, not been solved in 
going from Eq.(^) to Eq. fl2"6]) . However, it is possible to find a substitution for the 
spin wave function Eq.(pl[) that relates the calculation of Z[h] to an exactly solvable 



Brownian motion model. This method was already used by us in a previous paper fll 
where we considered the Hamiltonian Eq.([l7|) without interchain interactions (B = 0), 
in the context of disorder-induced solitons in a single polymer chain. To keep the present 
paper reasonably self-contained, we will briefly repeat the essential steps of this method, 
allowing now also for finite B values. Even though from a technical point of view, having 
B 7^ is not a major complication, the physical consequences are important, as it gives 
rise to the existence of LRO at finite temperature and (or) finite random magnetic field. 

We start by introducing the two time-dependent functions u{t) and v(t), such that 
the spin wave function Eq.Q2^)) reads: 

/ 

cosh( 



/ ___wu(t)\ \ 



tt(t) = u{t) 



(27) 



sinh(^) j 
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The initial condition, $?(t = 0) = | + ), requires that u(t = 0) = 1 and v (t = 0) = 0, and 
from Eq. (|25|) we obtain two coupled differential equations for u(t) and v(t), 

^ = u{t)e- pi cosh{v{t)), (28) 

and 

= -2(3 B - 2e~ (SI smh{v{t)) - 2f3h{t). (29) 

While Eg . (|28"D can be easily integrated, yielding 

u(t) = u{0) expje"^^ cosh(v (t')) dt'\ , (30) 

Eq.(p9|) can be interpreted as the Langevin equation for the velocity v of a particle 
undergoing a one-dimensional Brownian motion. The last term in Eq.(|29|) is the random 
force acting on the particle, while the first and the second terms describe, respectively, a 
constant (bias) force and a friction. 

In the absence of the random force (no random magnetic field), the velocity v ap- 
proaches a time- independent value for t — > oo: 

Woo = -arcsinh(/5 J Be /3/ ). (31) 

In the presence of the random force, v(t) becomes a stochastic variable with a time- 
dependent distribution, P(v,t), which may be obtained from the corresponding Fokker- 
Planck equation [18]. In the limit of long times (N — > oo), the time-independent equilib- 
rium distribution, P eq (v), is reached, obeying: 

-|(^W) + ?%M.^ . (32) 

Here, the first and second moments are given by, respectively, 
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ai(u) = -2(3B - 2e~ /3/ sinh(v(t)), (33) 

and 

a 2 = 4/? 2 / (h{t)h{t')) dt' = 4/? 2 e. (34) 
Jo 

The solution of the Fokker-Planck equation (j32|) is easily found to be: 

Peqiy) = Af exp | — z/t> — zcosh(t>)|, (35) 



where we introduced 



5 

0e' 



e -pi 



(3 2 e ' 

and the normalization coefficient M determined by: 



(36) 
(37) 



dvP eq {v) = 1. (38) 

The equilibrium distribution Eq.([35|) is centered around (Eq.([H])), and in the limit 
of vanishing fluctuations (e — > 0), it approaches a (^-distribution \P eq (v) — > 5(t> — «oo)]- 
Using our result for P eg (t> ), the average free energy of the RFIM can easily be obtained. 
We first note that Eq.(|30|) enables us to write the partition function Eq.(p6"D in terms of 
v{t) alone: 

Z[h] = 2^(N) = 2u(N) cosh = 

= 2 exp | e~ pi J* cosh(«(t)) dt j cosh (^f^J ■ (39) 

Omitting in the limit of long chains all terms that do not grow proportional to N, we 
can now replace the average of lnZ[/i] over the random field realizations by the average 



over P eq (v). This leads to 
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1 e _/3/ r+°° z K'(z) 

n < F > = " IT /_ cosh «" * = B ? ^ ■ (40) 

where K v (z) denotes the modified Bessel function of order v and K'„(z) its derivative 



with respect to z [p~9 ] . 



The analytical solution Eq.(KDJ) of the average free energy is the key result of this 
paper, from which all further results follow. In Fig. 1, this solution as a function of 
temperature is compared to results of a numerical simulation of the discrete RFIM at 
e = 0.01/ 2 for three values of the mean field B. Apart from the presence of the mean 



field, the simulation method is identical to the one described in Ref. Figure 1 shows 
that in the temperature region T <S I, the analytical solution is in excellent agreement 
with the simulation, except for very low temperatures, where the condition Eq.(|20"D is not 
met. It is observed that the analytical solution for the free energy reaches a maximum 
at T = T (e) ^ (indicated by the thin line). The meaning of T (e) will be discussed in 
Sec. Wn 



IV. ANALYSIS OF RESULTS 

Using Eq. flloD , it is straightforward to solve numerically the self-consistency equation 
(p|), from which the order parameter ((a)) is found through Eq.(^j). Figure 2 gives the 
thus obtained order parameter as a function of temperature and random-field strength for 
Ij_/I = 0.02. A low-temperature cut-off T (e) was used, as will be explained in detail in 



Sec. [IV B| . The thick curve in the ((ex)) = plane separates the region with LRO (inside 
this curve) from the one without LRO. From the solution for ((o - )), one obtains the critical 
temperature and the density of spin-flips through further numerical analysis. For certain 
parameter values v and z, however, our result Eq. ( fjO|) also allows for analytical solutions. 



13 



Here, two limiting cases should be distinguished. In Sec. [IV A| , we will consider the case 
of very weak disorder, where thermally induced spin-flips dominate the ones induced by 
the random magnetic field. In Sec. [IV B|, we consider the opposite limit. 



A. Dominant thermal spin-flips 

For a weak random magnetic field (e — > 0) the modified Bessel function in Eq.(|40|) 
has to be evaluated at large order and argument (u, z ^> 1). Using standard expansions 
T9fl , we obtain to first order in e: 



The first term in the right hand side of Eq. is the continuum approximation for the 
free energy of the discrete model Eq. (pi]) at e = 0. Under the conditions specified by 
Eqs.(|T9|) and fl20|) this term agrees with the exact expression, 

—F = — — In (cosh ([3 B) + ^(sinh ((3B)f + e" 2 ^) , (42) 



which can be obtained by diagonalization of the transfer-matrix Eq . ([14]) at hi = 0. 

Using Eq. ([41]) , the self-consistency equation ([□]) for the mean field B can be written 
in the form: 

Y x =V- ^e-^V , (43) 

where 

(f3B) 2 + e- 2/3/ ] _l . (44) 

At the transition temperature, T c (e), the order parameter vanishes. Therefore, by putting 
B = in Eq . (|43|) , we obtain an equation for T c (e): 
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In particular, at e = we obtain: 



exp 



2T C (0) " \T C (0) 



(46) 



This equation can be compared with the exact result for the critical temperature of the 



two-dimensional Ising model [20,21 



sinh ( — ^ sinh ( — ) = 1. 



For I± <C /, both Eq.(fl"6|) and Eq.pTj) give approximately: 



(47) 



Hi) 



(4f 



which shows that the chain mean-field approximation works well for the strongly 
anisotropic Ising model. 

From Eq. ([43]) we find that at weak disorder the transition temperature decreases 
linearly with e, 



T c (e) = T c (0) (1 - ae) 



(49) 



where a is: 



a 



I±(I + T c (0)) 



(50) 



Finally, the average density of spin-flips [see Eq.flP^D] for weak disorder is given by 

3-f 31 for T > T r (e) 



n, = < 



(51) 



. e ' m (ib + ep2y2 ) for T < Tc(e) 
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We note, that in the disordered phase (above T c (e)), no correction to the density of spin- 
flips linear in e occurs. This is related to the fact that in the weak disorder limit, the 
domain size R in Eq.([I]) is limited by the distance between thermally induced spin-flips, 
which is too small to allow for the creation of spin-flips by the random field. In this case 
the lowest-order correction is e 2 /? exp(/3I)/8. 

B. Dominant disorder-induced spin- flips 

We now turn to the limit exp [—/?/] <C f3 2 e, or equivalently: z < 1, where the thermal 
creation of spin-flips is negligible compared to creation by the random magnetic field. To 
lowest order in z, we then have in Eq. (j40|) : 



zK'jz) r(i + i/)(f) +r(i-i/ M2 



V 



v K u {z) T{l+v) (§) " - T(l 



V) 



(52) 



To find the critical temperature, at which i5 — > 0, it is sufficient to consider only small 
values of v oc B, where T(l ± v) ~ e Tlu (7 ~ 0.577 is Euler's constant), so that Eq.([52]) 
further reduces to: 

vK v {z) tanh(z/ln^)' 1 ' 

Using this result and the definitions Eq.(|36|) and Eq.(p7[) of v and z, respectively, the free 
energy Eq.fl40|) finally becomes: 

— (F) — . (54) 

N K ' I(T, e) tanh(x) 1 } 

Here, we introduced 

I(T,e) = I - 2rin(^r), (55) 

with 
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T = T (e) = V2e-7-2 e , (56) 

a characteristic temperature that depends on the strength of the random magnetic field. 
Furthermore, we defined x, 

X= ! JMl. (57) 
e 

At T = To, the free energy reaches a maximum (cf. Fig. 1), implying that the entropy 
vanishes. The general expression for the entropy in our model is: 

<s> = ~n~&f = ln • (58) 

which becomes negative for T < T . This unphysical behavior is also observed for isolated 
chains (Jj_ = 0) |lTfl . Its origin is the continuum approximation, which breaks down at 



low temperatures. As we discussed in Ref . [O , the problem arises due to the fact that in 



the continuum model a spin-flip can take any position on the chain and is not restricted 
to the lattice points of the original discrete model. The continuum model, therefore, gives 
incorrect results as soon as the thermal fluctuation of the spin-flip positions, Z(T) ~ T 2 /e, 
becomes less than one lattice constant. This happens at T < y/e, in agreement with 
Eq. (|56"D . The nature of T suggests that at this temperature the continuum model in 
fact describes the zero temperature discrete model. Using numerical simulations, we 
confirmed that for the one- dimensional RFIM this indeed is a reasonable identification 
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Using Eq. (|54|) , the self-consistency condition Eq.(|TT|) now takes the form: 

((a)) = ^ = ta^m>)"sinh 2 (x)- (59) 



For ((a)) 1, the non-trivial solution of Eq.(|59f) for the order parameter reads: 
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((a)) oc Jl{T,e) - i. (60) 



At the critical temperature T c (e), the left hand side of Eq. flBoD vanishes, yielding 

s( / - 2r -w h (^)) = 1 - (61) 

Figure 3 shows the critical temperature as a function of the random magnetic field 
strength for three values of the interchain interaction. The dashed curves are obtained by 
solving Eq. fl5"T|) . This yields meaningful results only for T (e) < T c (e) < 0.15J, where the 
upper limit arises from the fact that the zCl approximation made in deriving Eq.(|52]) 
starts to break down. The solid curves show the results obtained if we do not apply the 
small- z expansion, but instead numerically solve the self-consistency equation Eq . ([TT|) 
with the full free energy Eq . fl40|) . Clearly, over a large disorder interval, the approximate 
results give excellent agreement with the exact ones. Also presented (dots) are the results 
in the weak disorder limit (Eq.(f49|)). The dash-dotted curve in Fig. 3 indicates the lower 
temperature limit T (e), which we argued to correspond to T = in the discrete model. 

The critical disorder strength e c below which LRO exists in the ground state is a 
function of the interchain interaction and can be calculated from Eq.(^) by requiring 
that T c (e) = T (e). This yields 



^¥(i + fV^I, (^) 



where 



c = 2V2e-T" 2 ~ 0.78. (63) 

It is seen that e c is of the order of I±I. For a random field strength e <C e c , the ground 
state (T = To) of the anisotropic RFIM is nearly perfectly ordered, ((c)) ~ 1, as is clear 
from Fig. 2. 
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We now turn to the average density of spin-flips within the chains, which using 
Eqs.(|T2|) and (|54D, is obtained as: 



x 2 



ns ^ e) = iWJ)^wr (64) 

In the disordered phase (x = 0), Eq . (|64l) reduces to 

n s{T,e) = p ^ e y (65) 
In particular, the latter expression holds for the one- dimensional RFIM, which is dis- 



ordered at all temperatures Equation (^5]) was obtained previously in Ref. (TT 

where we discussed the density of kinks in isolated polymer chains. The temperature 
region where Eq.(|65|) is valid is specified by Eqs. fllTf ) and (p0[). However, as we showed in 



Ref. 0], the density of spin-flips at zero temperature in the discrete RFIM is very close 



to the continuum result Eq.(|65|) at T = T (e): 

n s (0,e) & 2 , (66) 

( I + cy/e) 

where c is defined by Eq. fl6"3|) . To lowest order in e, this result agrees with the estimate 
Eq.(^). The leading correction to Eq.(0) is of the order e 3//2 . 

For e < e c , LRO does exist and the zero temperature spin-flip density should be 
obtained by numerically solving the order parameter from Eq.(p9|) at T = T and substi- 
tuting the result into Eq.(|64|). For Ij_/I = 0.02 the resulting e-dependence of the spin-flip 
density is shown in Fig. 4. At this particular value of the interchain interaction, the crit- 
ical disorder strength is e c ~ 0.00 71 J 2 . It is observed from Fig. 4 that well below this 
critical value, the spin-flip density is strongly suppressed. In fact, for e e c , Eq.(|64j) 
reduces to 

n s (T , e) ~ — exp I I , (67) 
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which shows exponential suppression for small e. Equation fl6"T|) is represented in Fig. 4 
by the dashed curve. It should be noted that, strictly speaking, Eq.(^) is only valid 
for v = B/(f3e) <C 1 which was assumed in deriving Eq.fl53|) from Eq. fl52|) . However, 
straightforward analysis shows that for Jj_/J = 0.02 and T = T , this condition on v 
implies e ^> O(10~ 4 ), which holds for the overwhelming part of Fig. 4. 

V. DISCUSSION 

The results for the anisotropic Ising model obtained above admit a simple interpre- 
tation. To this end it is useful to distinguish between two kinds of spin-flips occuring in 
chains: kinks and antikinks. For kinks, the spins are positive to the left of the spin-flip 
and negative to the right. For antikinks, the opposite holds. Obviously, kinks and an- 
tikinks have the same creation energy, and along the chain always an alternation of kinks 
and antikinks occur. 

We first discuss the density of thermally induced spin-flips in the absence of a random 
magnetic field (Eq. (|51|) with e = 0). As we assumed the temperature to be much smaller 
than the spin-flip creation energy, J, the density of the thermally-induced spin-flips is 
small. In the disordered phase (above T c ), one can neglect the correlations between 
the positions of spin-flips in different chains and the density of the spin-flips (kinks and 
antikinks) is given by the Boltzmann formula: n s = e _/3/ . In the ordered phase, the 
interchain interactions bind kinks and antikinks within one chain into pairs, which is 
clear from the fact that n s oc e~ 2/SI below T c . As follows from Eq.(|3]) the binding energy 
V of such a pair is proportional to the distance R between kink and antikink: 

V(R) = I±R, (68) 

as over this distance the spin sign is opposite to ((ex)). Thus, below T c (0), the density of 
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the thermally induced kink-antikink pairs is: 

n pair ~ R{T)e~ m , (69) 

where the factor R(T) = T/(I±) - the average size of the pair at temperature T - comes 
from the sum over all possible distances between kinks and antikinks in the partition 
function. The density of spin- flips in the ordered phase equals 2n pa i r , in accordance with 

The transition between the ordered and disordered phases occurs when the kink- 
antikink pairs dissociate into free kinks and antikinks. Thus, at the critical temperature 
the average length of the kink-antikink pair is of the order of the average distance l/n s 
between spin-flips, 



R(T c (0)) ~ exp l^—j , (70) 



which explains the physical content of Eq. (|4q) [22 . 

Next we discuss the density of spin-flips in the disordered phase induced predom- 
inantly by the random magnetic field [see Eq. fl65|) 1. This equation coincides with the 
simple Imry and Ma estimate Eq.(§), except that the kink creation energy, J, is replaced 
by I(T,e) defined in Eq.([5^). This can be understood as follows: The spin-flips enter 
into the partition function with the weig ht: w(T) = Z(T)e _/3/ , where Z(T) ~ T 2 /e is the 
thermal fluctuation of the spin-flip position in the presence of the random magnetic field. 
We can now write the weight in the form: w(T) = e~^ I ^ T ' e \ where I(T, e) = I — Thil{T) 
is the spin-flip free energy, which also includes the entropy of the spin-flip location. It is 
clear that the free energy, rather than the "bare" kink creation energy, /, should enter 
into the correct expression Eq.(|65D for the density of spin-flips. 

Finally, we discuss the critical strength of the random magnetic field, given by Eq.fl62"D, 
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above which the system is disordered at all temperatures. In an isolated chain, spin-flips 
are induced by an arbitrarily weak random magnetic field [see Eq.(|65|)1, because the 
distances between the spin-flips can be made as large as is necessary to compensate for 
the spin-flip creation energy by the energy of the interaction with the random magnetic 
field. In the presence of interchain interactions, however, the distances between spin-flips 
are not allowed to be arbitrarily large, because of the potential Eq . (|68|) . The potential 
grows linearly with the distance R between kink and antikink, while the energy gain 
due to kink-antikink pair creation grows proportional to \/R. Thus, in the presence of 
interchain interactions, the equation describing the balance between the energy gain and 
energy loss due to the creation of a kink-antikink pair of size R [cf. Eq.(|l])], 



21 + I ± R ~ 2VeR , (71) 



only has a solution if e > e c ~ I±I, in accordance with Eq.(p2|). For e < e c the distances 
between disorder-induced spin-flips in isolated chains are very large, so that they are 
suppressed by the interchain interactions and the system at zero temperature is in the 
ordered state. If, on the other hand, e > e c the system is disordered at all temperatures. 

The dependence of the critical temperature on the disorder strength in the RFIM can 
be compared with that in another Ising-type model of disordered (spin-) Peierls systems 



considered in Ref. [22]. There it was assumed that impurities randomly cut the chains 
into finite segments, some of which, e.g., the ones with an odd number of units, contain 
at least one soliton. This was modeled by an ensemble of finite Ising chains whose length 
are taken from a Poisson distribution and which are subjected to a mean field and four 
different kinds of boundary conditions. For small impurity concentration x, the critical 
temperature in that model decreases linearly with x, just as in the RFIM [see Eq. ([49|) l. 
On the other hand, for large x, T c (x) ccl/x and at zero temperature the ground state is 
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ordered for all values of x. The absence of a critical concentration, above which LRO is 
destroyed, is related to the fact that in that model the positions of spin-flips (solitons) 
within the finite chains are not fixed. At T = the spin-flips are located near one of 
the chain ends to minimize the size of "wrong" domains. In the RFIM, however, the 
distances between neighboring spin-flips are governed by the energy balance Eq.(|7TJ) and 
cannot be made arbitrarily small, which results in the existence of a critical random field 
strength. 

VI. SUMMARY AND CONCLUSIONS 

In this paper, we studied the anisotropic RFIM, consisting of spin chains with in- 
terchain exchange interactions that are much weaker than the intrachain interactions. 
Treating the interchain interactions in the mean-field approximation, we found an ana- 
lytical solution Eq. ([40]) for the free energy. The crux of our approach lies in the fact that 
in this chain mean-field approximation the free energy of the continuum version of the 
model can be related to a one-dimensional Brownian motion. We found the stationary 
solution of the corresponding Fokker-Planck equation, which enabled us to perform the 
average of the free energy over the random magnetic field realizations. From this, we 
obtained the order parameter as a function of temperature and random field strength 
and, more specifically, the critical temperature below which LRO occurs as a function of 
the random field strength. We also calculated the density of spin-flips in the chains. As 
we decribed in Sec. 0, our results have a clear physical interpretation. 

The chain mean-field approximation restricts the validity of our results to d > 2, where 
thermodynamic quantities suffer much less from fluctuations than in lower dimensional 
systems. We believe that the critical temperature T c (e) given by Eq.(^TJ) describes the 
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phase transition in three dimensions qualitatively correct and is quantitatively correct 
for a strongly anisotropic realization of the RFIM. Clearly, our mean-field treatment of 
interchain interactions is not applicable in two dimensions, where LRO is lacking for all 
temperatures and (or) random-field strengths JTU . 



The one-dimensional RFIM had already been solved by us in the context of disorder- 
induced solitons in a single chain of the polymer trans-polyacetylene [|TTJ. Our present 
treatment of the anisotropic d- dimensional RFIM allows us to extend that study to 
include the soliton-antisoliton confinement that is imposed by electron hopping between 
neighboring polyacetylene chains in the bulk material. In particular, the exponential 
suppression of the density of spin-flips well below the critical random field strength (see 
Eq.flSTD) may be relevant to explain the absence of experimental signatures of these 



solitons. In a forthcoming paper we will focus more particularly on this application 
of the anisotropic RFIM and discuss various experimental observables. We finally note 
that if B is considered an external (rather than self-consistent) field, our model describes 
the lattice dimerization of a disordered polymer chain with a non-degenerate ground 
state ||12|| . B then is proportional to the energy difference between the two phases with 



different signs of the dimerization. 
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Figure captions 



Fig. 1. The average free energy of the anisotropic RFIM as a function of temperature 
for three values of the mean field: B = 0.005/, 0.01/, and 0.02/ (curves a, b, and 
c, respectively). Solid curves represent the analytical result Eq.(p3|) for the continuum 
approximation to the RFIM; dashed curves are obtained from numerical simulations of 
the RFIM. In all cases, the random field strength was taken e = 0.01/ 2 . The thin line at 
T = To gives the position of the maxima of the analytical solution. 

Fig. 2. The order parameter ((a)) of the ansitropic RFIM as a function of the random 
magnetic field strength e and temperature T (for T > T (e)). The plot was obtained by 
solving the self-consistency condition Eq . fllTD with the analytical solution Eq. (|40|) for the 
free energy. The interchain interaction was chosen such that I±/I = 0.02, in which case 
T c ~ 0.295/ (at e = 0) and e c = 0.0071/ 2 . At (T, e) = (0, 0), the order parameter reaches 
its maximum: ((a)) = 1. The thick curve in the ((cr)) = plane separates the (e, T) 
regions with and without LRO. 

Fig. 3. The critical temperature T c (e) of the anisotropic RFIM as a function of 
the random magnetic field strength for I±/I = 0.002, 0.01, and 0.02 (curves a, b, and c, 
respectively). Solid curves derive from solving Eq.(O) with the exact free energy Eq.([5U]). 
The dashed curves are obtained by solving Eq.(^l|), while the dots represent the small-e 
behaviour Eq.(||]). The dash-dotted curve shows the temperature T (e), below which the 
continuum approximation breaks down and which corresponds approximately to T = 
in the discrete RFIM. 
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Fig. 4. Density of spin-flips within the chains of the ansitropic RFIM as a function of 
the random magnetic field strength at T = Tq, i.e., in the ground state. The interchain 
interaction was chosen such that I±/I = 0.02. The dashed line shows the analytical result 
Eq.(§7p for e < e c ~ 0.0071/ 2 . 
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